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O ■ ABSTRACT 
^ ' 

Aims. To provide the observational constraints required to use the gravitationally lensed quasar SDSS J0924+0219 for the determination of Ho 
from the time delay method. We measure here the redshift of the lensing galaxy, we show the spectral variability of the source, and we resolve 
the lensed host galaxy of the source. 

rS I Methods. We present our VLT/FORSl deep spectroscopic observations of the lensed quasar SDSS J0924+0219, as well as archival 
HST/NICMOS and ACS images of the same object. The two-epoch spectra, obtained in the Multi Object Spectroscopy (MOS) mode, allow for 
very accurate flux calibration and spatial deconvolution. This strategy provides spectra for the lensing galaxy and for the quasar images A and 
B, free of any mutual light contamination. We deconvolve the HST images as well, which reveal a double Einstein ring. The mass distributions 
' in the lens, reconstructed in several ways, are compared. 
• • Results. We determine the redshift of the lensing galaxy in SDSS J0924+0219: zie„s= 0.394 ± 0.001. Only slight spectral variability is seen in 
. ^ continuum of quasar images A and B, while the C III], Mg II and Fe II emission lines display obvious changes. The flux ratio between 
, the quasar images A and B is the same in the emission lines and in the continuum. One of the Einstein rings found using deconvolution 
» I ■ corresponds to the lensed quasar host galaxy at z = 1.524 and a second bluer one, is the image either of a star-forming region in the host 
' galaxy, or of another unrelated lower redshift object. A broad range of lens models give a satisfactory fit to the data. However, they predict 
very different time delays, making SDSS J0924+0219 an object of particular interest for photometric monitoring. In addition, the lens models 
reconstructed using exclusively the constraints from the Einstein rings, or using exclusively the astrometry of the quasar images, are not 
compatible. This suggests that multipole-like structures play an important role in SDSS J0924-I-0219. 
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1. Introduction 

COSMOGRAIL is a multi-site optical monitoring campaign of 
lensed quasars. Following the original work by Refsdal ( 1964 1, 
its goal is to measure, with an accuracy close to one percent 
(Eigenbrod al. .2005< . the so-called time delay between the im- 
ages of most gravitationally lensed quasars. These time delays 



* Based on observations made with the ESO-VLT Unit Telescope 2 
Kueyen (Cerro Paranal, Chile; Program 074.A-0563, PI: G. Meylan) 
and on data obtained with the NASA/ESA Hubble Space Telescope 
(Program HST-GO-9744, PI: C. S. Kochanek) and extracted from the 
data archives at the Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. 



are used in combination with lens models and detailed observa- 
tions of individual systems to infer the value of the Hubble pa- 
rameter H(), independent of any standard candle (e.g., reviews 
by Courbin et al. 2002 Kochanek 2005a). 

The present work is devoted to the quadruply imaged 
quasar SDSS J0924+0219 (Inada et al. 2003) at z = 1.524, dis- 
covered in the course of the Sloan Digital Sky Survey (SDSS). 
This object is particularly interesting because of its anomalous 
image flux ratios, the origin of which is unclear. It has been 
argued that the faintest image of SDSS J0924-I-0219, which 
is located at a saddle point of the arrival-time surface, could 
be demagnified either from star microlensing (Schechter et al. 
i2004i Keeton et al. i2005 ) or subhalos microlensing (Kochanek 
& Dalal. l2004l . 
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Fig. 1. /?-band VLT image of SDSS J0924+0219, where objects are labeled following Inada et al. JMBt . The stars a, c, d, and 
e are used to compute the PSF spectrum (see text). Only stars a, d and e are used to derive the relative flux calibration between 
each MOS mask. The field of view is 3.4' x 3.4', North is to the top. East to the left. 



We analyse here our deep optical spectra of SDSS 
J0924+0219 obtained with the ESO Very Large Telescope 
(VLT). These spectra are used to: 1- measure the redshift of the 
lensing galaxy, 2- estimate the spectral variability of the quasar, 
3- measure the flux ratio between images A and B of SDSS 
J0924+0219, in the continuum and the broad emission Unes. 
Hubble Space Telescope (HST) ACS and NICMOS images 
from the STScI archives are deconvolved using the MCS al- 
gorithm (Magain et al. ll998> which unveils two Einstein rings. 
One of the rings corresponds to the host galaxy of the quasar 
source and is used to constrain the lens models. The second one 
is probably due to a star-forming region in the host galaxy of 
the quasar source or to another unrelated object. 

2. VLT Spectroscopy 

2.1. Observations 

Our spectroscopic observations of SDSS J0924H-0219 are part 
of a low dispersion spectroscopic survey aimed at measuring 
all unknown lens redshifts. They are acquired with the FOcal 



Reducer and low dispersion Spectrograph (FORSl), mounted 
on the ESO Very Large Telescope, used the MOS mode (Multi 
Object Spectroscopy) and the high resolution collimator. This 
configuration allows the simultaneous observation of a total of 
8 objects over a field of view of 3.4' x 3.4' with a pixel scale 
of 0.1" (Fig.Q. The G300V grism, used in combination with 
the GG435 order sorting filter, leads to the useful wavelength 
range 4450 A < A < 8650 A and to a scale of 2.69 A per pixel 
in the spectral direction. This setup has a spectral resolution 
R - A/AA ^ 300 at the central wavelength A - 5900 A, which 
translates in velocity space to Av = cAA/A ^ 1000 km s"'. 

The sHtlets of the MOS mask are all 19" long and 0.7" 
wide, which both avoids lateral contamination by the quasar 
image C and matches well the seeing values during the obser- 
vations. Four slits were centered on the foreground stars a, c, 
d, e, while a fifth slit is centered on images A and B of SDSS 
J0924-I-0219, after rotation of the mask to a suitable Position 
Angle (PA) (Fig. O. The spectra of the stars are used both to 
compute the reference Point Spread Function (PSF) needed for 
the deconvolution and to carry out a very accurate relative flux 
calibration. "Through-slit" images acquired just before expo- 
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Fig. 2. 7?-band images of SDSS J0924+0219. A short 30 sec 
exposure is shown on the left, where the quasar images A, B, C 
and D as well as the lensing galaxy, are indicated. The seeing 
is 0.37" in this image and the pixel scale is 0.10". The posi- 
tion of the 0.70" slitlets is also indicated. They correspond to 
three epochs of observations with very different seeings (see 
TableQ. The slit has not moved at all between the exposures, 
even when taken 15 days apart. 

Table 1. Journal of the VLT spectroscopic observations of 
SDSS J0924+0219. The seeing is measured on the spectrum 
of the PSF stars. The exposure time is 1400 s for each of the 6 
epochs. 



ID 


Date 


Seeing ["] 


Airmass 


Weather 


1 


14/01/2005 


0.66 


1.188 


Photometric 


2 


14/01/2005 


0.59 


1.150 


Photometric 


3 


14/01/2005 


0.46 


1.124 


Light clouds 


4 


01/02/2005 


0.83 


1.181 


Photometric 


5 


01/02/2005 


0.97 


1.146 


Light clouds 


6 


01/02/2005 


0.84 


1.126 


Light clouds 



sures # 1, # 3, # 4 in order to check the mask alignment are 
displayed in Fig. |2l 

2.2. Reduction and Deconvolution 

The spectra are bias subtracted and flat-fielded using IRAF^ 
The flat fields for each slitlet are created from 5 dome expo- 
sures, using cosmic ray rejection. They are normaUzed by av- 
eraging 60 lines along the spatial direction, rejecting the 20 
highest and 20 lowest pixels, then block replicating the result 
to match the physical size of the individual flat fields. 

Wavelength calibration is obtained from numerous emis- 
sion lines in the spectrum of Helium- Argon lamps. The wave- 
length solution is fitted in two dimension to each slitlet of 
the MOS mask. The fit uses a fifth-order Chebyshev polyno- 
mial along the spectral direction and a third-order Chebyshev 

' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



polynomial fit along the spatial direction. Each spectrum is in- 
terpolated following this fit, using a cubic interpolation. This 
procedure ensures that the sky lines are well aligned with the 
columns of the CCD after wavelength calibration. The wave- 
length solution with respect to the reference lines is found to 
be very good, with an rms scatter better than 0.2 A for all spec- 
tra. 

The sky background is then removed by fitting a second- 
order Chebyshev polynomial in the spatial direction to the areas 
of the spectrum that are not illuminated by the object. 

Finally, we perform the cosmic ray removal as follows. 
First, we shift the spectra in order to align them spatially (this 
shift is only a few tenths of a pixel). Second, we create a com- 
bined spectrum for each object from the 6 exposures, removing 
the 2 lower and 2 higher pixels, after applying appropriate flux 
scaling. The combined spectrum obtained in that way is cos- 
mic ray cleaned and used as a reference template to clean the 
individual spectra. We always check that neither the variable 
seeing, nor the variability of the quasar causes artificial loss of 
data pixels. 

Even though the seeing on most spectra is good, the lensing 
galaxy is close enough to the brightest quasar images A and B 
to be aff'ected by significant contamination from the wings of 
the PSF. For this reason, the spectral version of MCS decon- 
volution algorithm (Magain et al. 1998. Courbin et al. I2000> 
is used in order to separate the spectrum of the lensing galaxy 
from the spectra of the quasar images. The MCS algorithm uses 
the spatial information contained in the spectrum of a reference 
PSF, which is obtained from the slitlets positioned on the four 
isolated stars a, c, d, and e (Fig.[0. The final normalized PSF 
is a combination of the four PSF spectra. The six individual 
spectra are deconvolved separately, extracted, flux calibrated 
as explained in Section E31 and combined. The spectrum of the 
lensing galaxy is extracted from the "extended channel" of the 
deconvolved data, while the spectra of the quasar images are 
extracted from the "point-source channel" (see Courbin et al. 
I2OOOI . 

2.3. Flux Calibration 

Our absolute flux calibration is based on the spectrum of the 
spectrophotometric standard star Feige 66 taken on the night 
of 2005 January 16. The response function of the grism is de- 
termined for this single epoch. It is cross calibrated using stars 
observed in each MOS mask in order to obtain a very accurate 
calibration across all epochs. The spectra of four stars are dis- 
played in Fig. |3j without any deconvolution and having used 
a 4" aperture for extraction. We find significant differences in 
flux between the six epochs, that need to be corrected for. The 
main causes for these diff'erences are variable seeing and vari- 
able extinction due to thin cirrus during some of the observa- 
tions (Table 0. The effect of mask misalignment is excluded, 
as can be seen from the image-through-slit of Fig.|2] 

Assuming that the intrinsic flux of the foreground stars has 
not varied between the six exposures, and taking the data # 1 
of Tabled as a reference, we derive the flux ratio between this 
reference epoch and the six other dates, for each star. These 
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Fig. 3. The spectra of the foreground stars. The index on the 
right of each spectrum indicates the exposure number, follow- 
ing Table [2 Flux differences are mainly due to the presence of 
light clouds on observation dates # 3, # 5 and # 6. 



curves, fitted with a third-order polynomial, are shown in Fig. 
13 The corrections computed in this way are found to be very 
stable across the mask: the curves obtained for two different 
stars only showed slight oscillations with an amplitude below 
2%. This is also the accuracy of the flux correction between 
different epochs. A mean correction curve is then computed for 
each epoch from all stars, except star c which is much fainter 
than the others, and is applied to the deconvolved spectra of the 
quasars and of the lensing galaxy. 

3. Extracted Spectra 

3.1. The Lensing Galaxy 

The six deconvolved spectra of flie lensing galaxy are extracted, 
combined, and smoothed with a 5 A box (2 pixels). Fig. |5l 
shows the final one-dimensional spectrum, where the Ca II H 
& K absorption lines are obvious, as well as the 4000 A break, 
the G-band typical for CH absorption, the Mg band, and the H/3, 
and Fe II absorption lines. These features yield a mean redshift 
of ziens= 0.394 + 0.001, where the 1-cr error is the standard de- 
viation between all the measurements on the individual lines, 
divided by the square root of the number of lines used. We do 
not consider the 4000 A break in these calculations. This spec- 
troscopic redshift falls very close to the photometric estimate 
of z 



0.4 by Inada et al. J2003> . and agrees with the spectro- 
scopic redshift of Ofek et al. J2005> . In addition, the absence of 
emission lines confirms a gas-poor early-type galaxy. No trace 
of the quasar broad emission lines is seen in the spectrum of the 
lensing galaxy, indicative of an accurate decomposition of the 



Fig. 4. Flux ratios between Date # 1 and the 5 others, along 
with the third-order polynomial fits. We use the ratios of the 
3 stars: a, d and e to determine the mean correction applied 
to the quasar Star c, which is much fainter than the others, 
is excluded from the final calibration. The (small) parts of the 
spectra with strong atmospheric absorption are masked during 
the polynomial fit. The peak-to-peak differences between the 
ratios computed using star a, d and e, differ by less than 2%. 



data into the extended (lens) and point source (quasar images) 
channels. 



3.2. The Quasar Images 

The mean spectra of quasar images A and B are shown in 
Fig.il smoothed with a 5 A box. The Al III], Si III], C III], 
[Ne IV] and Mg II broad emission lines are clearly identi- 
fied. A Gaussian fit to these 5 lines yield a mean redshift of 
1 .524 + 0.001 for image A and 1 .524 + 0.002 for the fainter im- 
age B. The standard deviation between the fits to the individual 
lines, divided by the square root of the number of lines used, is 
taken as the error bar These results are in excellent agreement 
with the values obtained by Inada et al. <2003> . as well as the 
redshift from the SDSS database, who both report z = 1 .524. 

3.3. Variability of the Quasar Images 

The spectra of quasar images A and B are shown in Fig. 
for 2005 January 14 and February 1. These are the mean of 
the three spectra obtained on each date, smoothed with a 5 A 
box. Although the continuum shows little variation (only B has 
fadded slightly between our two observing dates), there are ob- 
vious changes in the broad emission lines of each quasar im- 
age. In image A, the red wing of the Mg II emission line has 
brightened, as well as the C II] emission line, while in image 
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Fig. 5. Spectrum of the lensing galaxy in SDSS J0924+0219, 
as obtained by combining the data for the 6 epochs, i.e., a total 
integration time of 8400s. The template spectrum of an ellipti- 
cal galaxy at z=0.394 is also shown for comparison (Kinney et 
al. 1 19961 1. All main stellar absorption Unes are well identified. 
Prospects for a future determination of the galaxy's velocity 
dispersion are therefore excellent. 



B, the center of the C III] emission line has become double 
peaked and the Fe II feature redwards of Mg II has fadded. 
A zoom on these lines is shown in Fig. |8] The line variations 
are already visible before averaging the 3 individual spectra at 
a given date and in the not-so-blended quasar images of the 
raw un-deconvolved spectra. We can therefore safely rule out 
any deconvolution artefacts due to PSF variations in the MOS 
mask. In addition, the residual images after deconvolution (see 
Courbin et al. 2000 for more details) are particularly good, in- 
dicative of little or no PSF variations across the slitlet mask. 

3.4. Image flux ratio 

Keeton et al. ( I2005 > have recently observed that the flux ra- 
tio between the images of SDSS J0924-I-0219 is different in 
the continuum and in the broad emission lines. In their slitless 
HST/ACS observations, the flux ratio between A and B is 2.60 
in the emission lines, and about 3.5 in the continuum, i.e., the 
emission lines are 30% different from the continuum. 

We plot the flux ratio between quasar image B and A as a 
function of wavelength at a given date (top panels in Figs. |9] 
and llOt . This ratio is close to flat, with some small differences 
in the broad emission lines. 

We construct the spectrum aB+fi and adjust the parameters 
using a linear least squares fit so that it matches the spectrum 
of quasar A. The result is shown in the middle panels of Figs.|5] 
andllOl Almost no trace of the emission lines are seen in the dif- 



Fig. 6. Spectra of the quasar images A and B of SDSS 
J0924-I-0219, as extracted from the deconvolved data. These 
figure shows the mean of the 6 spectra taken for the 6 epochs, 
after the flux calibration described in Section lOl 



ference spectra in the bottom panels of the figure. Our spectra 
indicate no strong differential amplification of the continuum 
and broad emission Unes in the components A and B of SDSS 
J0924-I-0219, and the small residual seen in the emission lines 
in the bottom panels of Figs.|5]and^|are an order of magnitude 
smaller than reported in Keeton et al. (I2005> . 

In the 15 days separating the observations, a has changed 
by only 2%. For both dates the residuals of the fit are almost 
perfectly flat, indicating no continuum change. Only asymmet- 
ric changes in the emission lines are seen. 

Finally, the flat flux ratio between image A and B shows 
that there is no significant extinction by interstellar dust in the 
lensing galaxy. 

3.5. Intrinsic variations vs. microlensing 

It is hard, with only two observing points, to infer the ori- 
gin of the spectral variations observed in SDSS J0924-1-0219. 
Nevertheless, we see rapid (15 days) and asymmetric changes 
in the emission lines of the quasar images, and no strong 
changes in the continuum. Intrinsic variations of quasars are 
usually stronger in the continuum than in the emission lines, 
and they are also longer than the two-week span we observe 
here. Such rapid variations due to microlensing have been seen 
in at least one other lensed quasar: HE 1104-1805 (Schechter 
et al.l2M3 l. SDSS J0924+0219 might be a second such case. 
Microlensing variability is supported by the photometric broad- 
band data by Kochanek et al. ( I2005bl) . showing that A and B 
have very different light curves that are hard to match even after 
shifting them by the expected time delay. However, microlens- 



6 



A. Eigenbrod et al.: COSMOGRAIL II: SDSS J0924+0219: lens redshift, quasar source spectra, and Einstein rings 




Fig. 7. The spectra of images A and B on 14 January and 1 
February 2005 show a stable continuum for both images, but 
the broad emission lines do vary on a time-scale of two weeks 
(see Fig. 



ing usually acts on the continuum rather than on the emission 
lines of quasar spectra, because of the much smaller size of 
the continuum region. Differential amplification of the contin- 
uum relative to the emission lines, as observed by Keeton et 
al. (I2005> . would be a strong support to the microlensing hy- 
pothesis. Our spectra do not show such a differential amplifi- 
cation, but we note that our wavelength range is very different 
from that of Keeton et al. J2005t and that they observed in May 
2005, i.e., 3 months after our observations. 

Assuming microlensing is the correct interpretation of the 
data, its strength depends upon the scale-size of the source, 
with smaller sources being more susceptible to large magni- 
fication (e.g. Wambsganss & Paczvnski .1991) . The continuum 
emitting region and the broad-line region (BLR) of a quasar 
can appear small enough to undergo significant magnifications. 
The limiting source size for microlensing to occur is given by 
the Einstein radius projected onto the source plane. This means 
that only structures in the source with sizes comparable to or 
smaller than this radius will experience appreciable amplifica- 
tion. The Einstein radius, projected onto the source plane for 
microlenses with masses in the range 0.1 Mq < M < lOM© 
is 7 < < 70 light-days for a cosmology with Q„j = 0.3, 
Qa = 0.7 andhioo=0.65. 

Kaspi et al. (20001 derived sizes for active galaxy nuclei 
from reverberation mapping of the Balmer lines. As a func- 
tion of intrinsic luminosity, they found a global scaling of the 
broad-line region (BLR) ranging from approximately 1 to 300 
light days, which compares well with the Einstein radius of the 
microlenses in the lensing galaxy of SDSS J0924+0219. 



Fig. 8. Enlargements of Fig. comparing the broad emission 
lines of images A and B on 14 January (sohd curve) and 1 
February 2005 (dotted curve). Obvious variations are seen in 
the red wing of the Mg II in image A, in the center of the C III] 
in image B. The Fe II feature redwards of Mg II in image B has 
also changed by 20%. These variations are asymmetric about 
the center of the lines. The asymmetry is different in C III] and 
Mg II. 



The observations also reveal that the broad emission lines 
and the continuum do not vary on the same time scale. Indeed, 
the continuum of image A remains constant over the 15 -day 
time span of the observations, while the broad emission lines 
vary. 

Detailed microlensing simulations by Lewis & Ibata J2004> 
show that the correlation between the magnification of the BLR 
and the continuum source exists, but is weak. Hence variations 
in the broad emission lines need not be accompanied by vari- 
ations in the continuum. This argument has been confirmed 
through observations of other gravitationally lensed quasars 
(Chaitas et al. 2002, Richai'ds et al. 120041) . 

Another observational fact that needs some enlightening is 
the asymmetric amplification of the broad emission lines (see 
Fig. Such an amplification occurs for the C II] and Mg II 
emission lines in the spectra of image A. The red wings of these 
lines are significantly more amplified than the blue ones. An ex- 
planation for this is given by Abajas et al. (2002 1 and Lewis & 
Ibata ( 2004i) . who show that emission lines can be affected by 
substantial centroid shifts and modification of the line profile. 
Asymmetric modification of the line profile can be indicative 
of a rotating source. Microlensing of the part of the BLR that is 
rotating away from us would then explain the observed asym- 
metric line amplifications. This would imply that a microlens- 
ing caustic is passing at the edge of the broad line region, and 
is far enough from the continuum to leave it unaffected. 
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Fig. 9. Comparison between the spectra of images A and B 
taken on 14 January 2005. The top panel shows the dimen- 
sionless ratio B/A. The mean ratio is 0.32. In the middle panel, 
a first-order polynomial aB+/3 is fit to the spectra of image A. 
The best fit is obtained with a = 2.80 + 0.05 and /3 = 0.37. The 
difference in flux between A and the fitted aB+/3 polynomial 
is displayed in the bottom panel, and does not exceed a few 
percent of the flux. 



4. HST Imaging 

Optical and near-IR images of SDSS J0924+0219 are available 
from the HST archive in the F555W, F814W and F160W fil- 
ters. The F555W and F814W observations have been obtained 
on 18 November 2003 with the Advanced Camera for Surveys 
(ACS) and the Wide Field Channel (WFC). The F555W data 
consist of two dithered 1094 s exposures, each one being split 
in two (CRSPLIT=2) in order to remove cosmic rays. Two con- 
secutive 1 148 s exposures have been taken through the F814W 
filter, one hour later, again splitting the exposure time in two. 
Finally, the NICMOS2 observations, taken on 2003 November 
23, consist of 8 dithered exposures, for a total of 5312 s. The 
5-day period separating the optical and near-IR observations is 
of the order of the expected time delay between images A and 
B of the quasar. 

4.1. Image Deconvolution 

The MCS algorithm (Magain et al. ll998l l is used to deconvolve 
all images. This algorithm sharpens the images and preserves 
the flux of the original data. It also decomposes the data into a 
set of analytical point sources (the quasar images) and a numer- 
ical "extended channel" which contains all the features other 
than point sources, i.e., the lensing galaxy and the Einstein 
ring. All images are rebinned to a common pixel scale prior 
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Fig. 10. Same as in Fig. |9] but for the spectra taken on 1 
February 2005. The mean B/A ratio is 0.31, and the best fit 
of image A is obtained with a - 2.86 + 0.05 and /? = 0.43. 



to deconvolution and combined with cosmic ray rejection. The 
reference image adopted to carry out the whole deconvolution 
work is the first image taken through the F814W filter, i.e., im- 
age j8oi33031 in the HST archive. The position angle of this 
reference image relative to the North is PA=-2.67°. All the as- 
trometry in the following is given in the coordinate system of 
this image. The data used here are the pipeline-drizzled images 
available from the archive. The pixel scale in the deconvolved 
image is half that of the original image, i.e., 0.025"x0.025". 
The spatial resolution is the same in all deconvolved images, 
i.e., 0.05" Full- Width-Half-Maximum (FWHM). 

As the HST PSF has significant spatial variations across the 
field of view, stars located far away from SDSS J0924-I-0219 on 
the plane of the sky are not ideal for use in the image decon- 
volution. To circumvent this problem we have devised an iter- 
ative procedure. We first deconvolve the images with a fixed 
PSF, directly measured from stars. This gives a deconvolved 
image of the lens and Einstein ring, that we reconvolve with 
the PSF and subtract from the original data. A second PSF 
is re-computed from this new lens- and ring-subtracted im- 
age, directly from the quasar images, following the procedure 
described in Magain et al. (2005 1. This is similar to a blind- 
deconvolution, where the PSF is modified during the deconvo- 
lution process. A new deconvolved image is created with the 
improved PSF, as well as a new lens- and ring-subtracted im- 
age. We repeat 4 times in a row the procedure until the residual 
map (Magain et al. 1998 Courbin et al 1998 1 is flat and in aver- 
age equal to 1 cr after deconvolution, i.e., until the deconvolved 
image becomes compatible with the data in the;^f^ sense. 
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Fig. 11. Left: composite HST image using the observations through the F555W, F814W and F160W fihers. The resolution is 
respectively 0.10" in F555W and F814W, and 0.15" in F160W. Right: deconvolved image. It has a pixel size of 0.025" and a 
resolution of 0.05". The lensed host galaxy of the quasar is clearly seen as red arcs well centered on the quasar images. A second 
set of bluer arcs inside and outside the area delimited by the red arcs is also revealed. The field of view is 3.0" on a side. The 
image is slightly rotated relative to North, which is at PA=-2.67°. East is to the left. The white square shows the position of the 
perturber found for the SIE and NFW models of Section l5T2l 



Table 2. Astrometry of SDSS J0924+0219 and flux ratio be- 
tween the images. All positions are given relative to the lensing 
galaxy in the coordinate system of our reference HST image 
j8oi33031. The 1-cr error bar on the astrometry is 0.005", 
mainly dominated by the error on the position of the lensing 
galaxy. The error bar on the flux ratio is of the order of 10% for 
images B, C and 20% for image D, and includes the systematic 
errors due to the presence of the Einstein ring (see text). 



Object 


X 


Y 


F555W 


F814W 


F160W 




(") 


(") 








Lens 


+0.000 


+0.000 








A 


-0.185 


+0.859 


1.00 


1.00 


1.00 


B 


-0.246 


-0.948 


0.51 


0.46 


0.44 


C 


+0.782 


+0.178 


0.39 


0.34 


0.32 


D 


-0.727 


+0.430 


0.06 


0.06 


0.03 



4.2. Results 

The deconvolved images through the three filters are shown 
in Fig. as a colour composite image. Two sets of arcs are 
clearly seen, corresponding to the host galaxy of the source 
quasar, and to a bluer object not centered on the images of the 
quasar. This arc is well explained by a second lensed source 
(see Section l5T2t which is either a star-forming region in the 
source, or another unrelated object. 



Instead of using the conventional version of the MCS de- 
convolution algorithm, we use a version that involves a semi- 
analytical model for the lensing galaxy. In this procedure, the 
analytical component of the lensing galaxy is either a two- 
dimensional exponential disk, or a de Vaucouleurs profile. All 
slight departures from these two profiles are modeled in the 
form of a numerical array of pixels which includes the arcs as 
well. 

In all bands, we find that an exponential disk fits the data 
much better than a de Vaucouleurs profile, which is surpris- 
ing for an elliptical galaxy, as indicated by the VLT spectra. 
Table |2 gives a summary of our astrometry, relative to the 
center of the fitted exponential disk. The mean position an- 
gle of the lensing galaxy, in the orientation of the HST im- 
age, is PA - -61.3 + 0.5° (positive angles relative to the 
North in a counter-clockwise sense) and the mean ellipticity 
is e = 0.12 + 0.01, where the error bars are the dispersions 
between the measurements in the three filters. We define the 
ellipticity as e = 1 - b/a, where a and b are the semi-major and 
semi-minor axis respectively. Note that although the formal er- 
ror on the lens ellipticity and PA is small, the data show evi- 
dence for isophote twisting. The eff'ective radius of the galaxy 
is R, = 0.50 + 0.05". 

The flux ratios of the quasar images are derived from the 
deconvolved images. The MCS algorithm provides the user 
with the intensities of all point sources in the image, decon- 
taminated from the light of the extended features, such as the 
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SIE de Vaucouleurs NFW 




Fig. 12. The three plots give the reduced x'' as a function of lens ellipticity e and external shear j for the three analytic models 
used in the LENSMODEL package. No constraint is used on the image flux ratios. The contours correspond to the 1, 2 and 3-cr, 
confidence levels. The degeneracy between ellipticity and shear is clear. Only the NFW models are (marginally) compatible with 
no external shear at all, as also suggested by the semi-linear inversion of Section l5^ The black square in each panel indicated 
the best fit model, which parameters are given in Table|3l 



Table 3. Best-fit parametric models for SDSS J0924-I-0219, obtained with the LENSMODEL package (Keeton l200Tt . The po- 
sition angles of the lens and of the external shear Qy are given in degrees, positive angles being counted counter-clockwise 
relative to the North. The coordinates (x, y) of the centres of the models are given in arcseconds, and the time delays Af are ex- 
pressed in days relative to the leading image B. The extreme values for the time delays within the smallest 1-cr region of Fig. 1121 
are also given. We adopt a (0,„, Qa) = (0.3, 0.7) cosmology and hi()o=0.65. All models have one degree of freedom. 



Model Parameters y) e 6c y By At^B Afcs Atos 



SIE 


V 


= 0.87 


(- 


-0.003, 0.002) 


0.13 


-73.1 


0.053 


65.4 




9 110.4 


6 2" 


0.91 


de Vaucouleurs 


b -- 

Re 


= 2.64 
= 0.50 


(- 


-0.004, 0.002) 


0.16 


-70.1 


0.096 


77.3 


8.6-J 


1 Q Q 14.4 


949.7 


1.41 


NFW 


rs 


= 0.70 
= 1.10 


(- 


-0.003, 0.001) 


0.10 


-72.0 


0.047 


65.4 


4 qS O 
^•^3.6 


7.8i- 


5 4" 


0.72 



ring in SDSS J0924-I-0219 and the lensing galaxy. The error 
on the quasar flux ratio is dominated by the contamination by 
the Einstein ring. If the intensity of a quasar image is overes- 
timated, this will create a "hole" in the deconvolved Einstein 
ring at the quasar image position. If it is underestimated, the 
local at the position of the quasar image will become much 
larger than 1 cr. The flux ratios in Table|2lare taken as the ones 
giving at the same time a continuous Einstein ring without any 
"hole", and leading to a good^^, close to 1, at the position of 
the quasar images. The error bars quoted in Tableware taken 
as the difference between these two extreme solutions, divided 
by 2. They include both the random and systematic errors. 

5. Modeling 

Constraining the mass distribution in SDSS J0924-I-0219 is 
not trivial. Firstly, we do not have access to the true image 
magnifications due to contamination by microlensing and sec- 



ondly, the light distribution of the lensing galaxy is not very 
well constrained. The ellipticity and position angle of the lens 
change with surface brightness, indicative of isophote twisting. 
Measuring the faintest isophotes on the HST data leads to PA 
- -25°, as is adopted by Keeton et al. J2005> in his models. 
However, brighter isophotes and fitting of a PSF-deconvolved 
exponential disk profile yields PA - -61.3°. 

As a blind test for the shape of the mass distribution un- 
derlying the light distribution, and without using any constraint 
on the ellipticity or PA of the lens, we use the non-parametric 
models of Saha & Williams (2004i. Fitting only the image po- 
sitions returns an asymmetric lens whose major axis is aligned 
approximately East- West (i.e., PA - 90°). Given the discrep- 
ancy between this simple model and the observed light distri- 
bution, we test in the following a range of models with differing 
levels of observational constraints, in order to predict time de- 
lays. 
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X (arcsec) X (arcsec) 

Fig. 13. Annular mask applied to the F160W (left) and F555W 
(right) data with point sources masked out. The annulus in the 
F555W image is shifted by 0.1" to the left and 0.2" to the top 
with respect to the F160W image, to properly encompass the 
blue arc seen in Fig.fTTI 




0.5 1 1,5 2 0,5 1 1,5 2 



X (orcsec) X (orcsec) 

Fig. 14. Reconstructed source from F160W data (top left) and 
its lensed image (bottom left). A second source lying on the 
rightmost cusp caustic (top right) is reconstructed from the 
F555W image corresponding to the blue arc (bottom right). 



5.1. Parametric Models 
5.1 .1 . Using the flux ratios 

The LENSMODEL package (Keeton jSOOTli is used to carry 
out an analytical modeling of the lensing galaxy. Three lens- 
ing galaxy models are considered: the Singular Isothermal 
Elhpsoid (SIE), the Navarro, Frenk & White profile 
(NFW), and the de Vaucouleurs J1948> profile. In a first at- 
tempt, we constrain these models with the lensing galaxy po- 
sition, the relative positions of the lensed images (Tabled) and 
their flux ratios (taken as the mean of the ratios measured in the 



three F555W, F814W, F160W filters). If no external shear is in- 
cluded in the models, we find a lens ellipticity of e ^ 0.3 with 
a P.A. 6e ^ 85° and an associated ^ 200. The ellipticity and 
PA agree well with the models obtained from the semi-linear 
inversion method of Warren & Dye J2003> (see Section]^}. 

Next, we include external shear to the model. The lens po- 
sition angle 6^, coordinates, and ellipticity agree better with the 
measured values in the HST images. The;^^^ values remain bad 
(X^ - 30), although improved with respect to the models with- 
out external shear. The shear orientation is 6y ~ 60° which is 
about in the direction of a bright galaxy located 9.5" away from 
SDSS J0924-I-0219 and at PA = 53°. 

The main contribution to the total is the anomalous flux 
ratios between the images of SDSS J0924-I-0219. In particular, 
the extreme flux ratio between image A and image D of ~ 15, 
when these two images are predicted to have approximately the 
same brightness. This is not surprising because of the evidence 
of microlensing in image A (Sect. 13. 5> and of possible milli- 
lensing induced by massive substructures. This lead us to the 
considerations presented in the next section. 



5.1 .2. Discarding tlie flux ratios 

The modeling is similar to that of Sect. 15.1.11 External shear 
is included but the flux ratios are discarded. In order to use 
only models that have one degrees of freedom (DOE), we have 
fixed the effective radius of the de Vaucouleurs model to the 
observed value. Given the number of observational constraints, 
the NEW model would have zero DOE if all its parameters were 
left free during the fit. We have therefore fixed the orientation 
of the external shear in this model to the value we found in the 
SIE+shear model. The best fit models are presented in Table|3l 
with (reduced) improved to values close to 1. 

We map the lens ellipticity vs. external shear plane in or- 
der to estimate the degree of degeneracy between these two 
parameters. The results are displayed in Eig. It is imme- 
diatly seen that the 1-cr ellipses of the different models only 
marginally overlap. This is confirmed by the time delay values 
summarized in Table |3] where we also give the extreme val- 
ues of the time delays within the 68% confidence interval. The 
minimum difference between the extreme time delays predicted 
with a constant mass-to-light ratio galaxy (de Vaucouleurs) and 
by the more physically plausible SIE model is about 8%. Since 
the error measurement on the time delay propagates linearly 
in the error budget, even a rough estimate of the three time 
delays in SDSS J0924H-0219, with 8% accuracy will already 
allow to discriminate efficiently between flat M/L models and 
SIE. Distinguishing between SIE and NEW is more difficult as 
time delays predicted by NEW models differ by only 1 % from 
the SIE time delays. Such an accuracy will be hard to reach in 
SDSS J0924-I-0219, that has short time delays and a short vis- 
ibility period given its equatorial position on the plane of the 
sky (see Eigenbrod et al. 2005). 
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Table 4. Minimized lens model parameters and corresponding Model parameters are: a-q = mass normalization in arbitrary 
units, (x, y) = offset of lens model centre from lens optical axis in arcseconds, e - ellipticity, y = external shear, 6g and 6y 
= PA in degrees counted counter-clockwise from North. In the case of the NFW, the scale radius is held fixed at 6" in the 
minimization. The third column gives the number of degrees of freedom (NDOF). Subscript 'i>' refers to the secondary SIS in 
the dual component models (see text). 



Model 


/V tnm 


NDOF 


Minimized parameters 






SIE 


4280 


3975 


Ko = 100.0, (x,y) = (0.02,0.04), e 


= 0.270, e, = 86.0 




NFW 


4011 


3974 


Ko = 100.0, (x,y) = (0.06,0.06), e 


= 0.187, 0, = 84.9 




Dual SIS 


4385 


3974 


Ko = 49.2, (x,y) = (0.00, 0.28), kqi, 


= 51.6, (x,y)h = (-0.06,-0.33) 




SIE+SIS 


4247 


3972 


Ko = 99.4, (x,y) = (0.04, 0.04), e = 


0.265, 0, = 85.1, KOI, = 2.1, (x,y)i, 


= (-0.79,-0.03) 


NFW+SIS 


3971 


3971 


Ko = 98.0, (x,y) = (0.05, 0.08), e = 


0.206, e, = 83.1, Kob = 2.8, (x,y% 


= (-0.80,-0.09) 


NFW+7 


3992 


3972 


Ko = 100.0, (x,y) = (0.06,0.06), e 


= 0.168, e, = 86.0, y = 0.010, Oy = 


78.3 



5.2. Using the Arcs of the Lensed Sources 

The HST images of SDSS J0924+0219 reveal two sets of arcs. 
One is prominent in the near-IR (in red in Fig. II It and is well 
centered on the quasar images. It is the lensed image of the 
quasar host galaxy. A second set of bluer arcs is best seen in 
F555W. It is off-centered with respect to the quasar images, 
indicating either a companion to the quasar host, or an inde- 
pendent intervening object along the line of sight. 

We apply the semi-linear inversion method of Warren & 
Dye ( 2003 1 to the arcs observed in the F555W and F160W data. 
The method incorporates a linear matrix inversion to obtain the 
source surface brightness distribution that gives the best fit to 
the observed lensed image for a given lens model. This linear 
step is carried out per trial lens parametrisation in a standard 
non-linear search for the global best fit. 

Dye & Warren (I2005t successfully apply this technique to 
the Einstein ring system 0047-2808. They demonstrate that the 
extra constraints provided by the image of the ring results in 
smaller errors on the reconstructed lens model, compared to 
using only the centroids of the principal images as constraints 
in this system. 

In the case of 0047-2808, the source is a star form- 
ing galaxy without any point-like emission whereas the im- 
age of SDSS J0924-H0219 is clearly dominated by the QSO's 
central point source. To prevent the reconstruction of SDSS 
J0924H-0219 from being dominated by the point source and be- 
cause in this section only the reconstruction of the QSO host 
emission is of interest, we masked out the four point source 
images in the F555W and F160W data supplied to the semi- 
linear inversion code. The astrometry of the quasar images is 
not used as a constraint. Fig.ll3lshows the masked ring images. 

5.2.1. Reconstruction results 

The deconvolved F160W and F555W data are reconstructed 
with 6 different parametric lens models. Three of these are sin- 
gle mass component models: the singular isothermal ellipsoid 
(SIE), the elliptical NFW, and the elliptical NFW with external 
shear The remaining three test for asymmetry in the lens model 
by including a secondary singular isothermal sphere (SIS) mass 
component that is also free to move around in the lens plane 
and vary in normalization in the minimization. These models 



are the dual SIS model, the SIE-nSIS model and the NFW+SIS 
model. 

Since the F160W data have the highest signal to noise arcs, 
we base our lens modeling on these data and applied our overall 
best fit model to the F555W data to reconstruct the source. In 
all cases, we reconstruct with a 0.5" x 0.5" source plane com- 
prising 10 X 10 pixels. The reconstruction is not regularised, 
except in Fig. O where first order regularisation (see Warren 
& Dve l2003t is applied to enhance visualization of the source. 

Table|3]lists the minimized parameters for each model and 
the coiTesponding values ofx"^- The SIE-nSIS and NFW+SIS 
models clearly fare better than their single component coun- 
terparts, implying the lens is asymmetric. For the SIE+SIS, 
a decrease in of Ax^ = 33 for 3 fewer degrees of free- 
dom has a significance of 5.1 cr. The decrease of Ax^ - 40 
for the NFW+SIS has a significance of 5.7 cr. Both models 
consistently place the secondary SIS mass component around 
(-0.80", -0.05") with a normalization of only ~ 2.5% of the 
main component. 

Interestingly, the elliptical models listed in Table 01 have 
ellipticities close to those obtained with the LENSMODEL 
software, when no external shear is considered. When exter- 
nal shear is added to the NFW model, we do indeed obtain a 
significantly better fit compared to the NFW on its own, but 
the results differ from those listed in Tabled While the elUp- 
ticity remains almost the same as in Tabled its PA differs by 
approximately 25°. Moreover, we find a ten times smaller am- 
plitude for the shear using the semi-linear inversion than using 
LENSMODEL. Note, however, that the observed quasar image 
astrometry is used in the LENSMODEL analysis, whereas it 
is not in the present semi-linear inversion. If we use the lens 
model found by the semi-linear inversion to predict the posi- 
tion of the quasar images, we find poor agreement between 
the predicted and the measured positions. The global, large 
scale shape of the lens found by the semi-linear inversion is 
well adapted to model the Einstein rings, which are very sensi- 
tive to azimuthal asymmetry in the lens, but additional smaller 
scale structures are needed to slightly modify the positions of 
the quasar images and make them compatible with the mea- 
sured astrometry. The disagreement between the astrometry 
predicted by LENSMODEL and the one predicted by the semi- 
linear inversion adds support to the presence of multipole-type 
substructures in the lens (e.g., Congdon & Keeton 2005J . 
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The top left plot in Fig.^^shows the reconstructed source 
corresponding to the best fit NFW+SIS model for the F160W 
data. The observed arcs are explained by a single QSO host 
galaxy. Note that in this figure, purely to aid visualization, we 
have regularised the solution and plotted the surface brightness 
with a pixel scale half that used in the quantitative reconstruc- 
tion. The bottom left corner of Fig.^^shows the image of the 
reconstructed source lensed by the best fit NFW+SIS model. 

We then take the best fit NFW+SIS model in order to recon- 
struct the F555W data shown on the right in Fig.^] Note that 
the annular mask is shifted slightly compared to the F160W 
data, to properly encompass the blue arc. The reconstructed 
source and corresponding lensed image are shown on the right 
hand side of Fig. 1141 

There are two distinct sources now visible. The QSO host 
identified previously has again been reconstructed. This is be- 
cause its dominant image, the bright arc in the top left quadrant 
of the ring, is still present in the F555W data. A second source, 
more diffiise and lying on the rightmost cusp caustic is also 
visible. This second source is responsible for the blue arcs. 

The redshift of the second source remains unknown. It 
could be a star forming object/region lying 0.2" -D, - 1200 /i ^5,) 
pc away from the quasar, i.e., it would be part of the host 
galaxy. It is, however, not excluded that this second source is 
at a different redshift than the quasar, e.g. located between the 
quasar and the lens, as it is bluer than the quasar host galaxy. 
If the latter is true, SDSS J0924+0219 might be a unique ob- 
ject to break the mass sheet degeneracy. Unfortunately, the lens 
modeling alone, does not allow to infer a redshift estimate. 

5.3. Note on the different types of models 

The two methods used above differ in several respects. 
LENSMODEL has a limited number of free parameters but 
uses only the constraints on the astrometry of the quasar im- 
ages. While a qualitative representation of the lensed host 
galaxy of the quasar source can be attempted, the method does 
not allow a genuine fitting of the Einstein rings assuming a 
(simplified) shape for the quasar host. 

The semi-linear inversion carries out a direct reconstruction 
of the lensed source as a whole, where each pixel of the HST 
image is a free parameter. As the quasar images largely domi- 
nate the total flux of the source, they need to be masked before 
the reconstruction. For this reason it is not possible with this 
method, at the present stage of its development, to constrain the 
lens model using simultaneously the astrometry of the quasar 
images and the detailed shape of the Einstein rings. 

Although the two methods used in the present work are fun- 
damentally different and although they use very different obser- 
vational constraints, they agree on the necessity to bring extra 
mass near image D of SDSS J0924+0219. Smooth lenses like 
the ones implemented in LENSMODEL have PAs that differ by 
10° from the one measured in the HST images. In the orienta- 
tion of Fig. ^ the mass distribution found by LENSMODEL 
is closer to horizontal (PA= -90°) than the light distribution, 
hence giving larger masses next to image D. In the semi hnear 



inversion, the optimal position found for the SIS perturber is 
also close to image D. 

Given the above discussion, the poor determination of the 
lens PA is a main limitation to the interpretation of the time de- 
lays in SDSS J0924+0219. An alternative route is to determine 
the dynamical rotation axis of the lens, a challenge which is 
now within the reach of integral field spectroscopy with large 
telescopes and adaptive optics. 

6. Conclusions 

We have spatially deconvolved deep sharp VLT/FORS 1 MOS 
spectra of SDSS J0924+0219, and measured the redshift of the 
lensing galaxy, ziens= 0.394 ± 0.001, from numerous stellar ab- 
sorption lines. The spectrum beautifully matches the elliptical 
galaxy template of Kinney et al. ( 1996 1. 

The flux ratio between image A and B is Fa/Fb - 2.80 + 
0.05 on 2005 January 14, and Fa/Fb = 2.86 ± 0.05 on 2005 
February 1, i.e., it has not changed between the two dates given 
the uncertainties on the flux ratios (Table For each date, 
this ratio is mostly the same in the continuum and in the broad 
emission lines of the quasar images A and B. This may seem 
in contradiction with Keeton et al. (2005 1 who see differential 
amplification of the continuum relative to the lines, but our ob- 
serving dates and setup are very different from theirs. 

While the continuum of images A and B has not changed in 
15 days, there are obvious and asymmetric changes in some of 
the quasar broad emission lines. Microlensing of both A and B 
is compatible with this, although somewhat ad hoc assumptions 
must be done on the position of the microcaustics relative to 
the quasar, as well as on the relative sizes of the continuum and 
broad line regions. 

Deep HST imaging reveals two sets of arcs. One corre- 
sponds to the red lensed host galaxy of the quasar and de- 
fined an Einstein ring connecting the quasar images. The other, 
fainter and bluer is off-centered with respect to the quasar im- 
ages. It is either a star-forming region in the quasar source host 
galaxy, or another intervening object. 

The lens ellipticity and PA measured in the HST images are 
hard to reconcile with simple models without external shear. 
The model fits improve when external shear is added, even 
though the predicted PA differs from the measured one by ap- 
proximately 25°. 

Models of Section l?!2l involving an additional small (SIS) 
structure to the main lens always place it along the East- West 
axis, about 0.8" to the East of the main lens, i.e., towards the 
demagnified image D. In addition, the models reconstructed us- 
ing only the Einstein rings do not predict the correct astrome- 
try for the quasar images. Einstein rings constrain the overall, 
large scale of the lens. Small deviations from this large scale 
shape are needed to match the quasar images astrometry. The 
discrepancy between the models using the rings and the ones 
using the quasar image positions therefore adds support to the 
presence of multipole-like substructures in the lens of SDSS 
J0924+0219. 

Finally, the range of time delays predicted by the different 
lens models is large and is very sensitive to the presence of 
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external shear and to the determination of the main lens ellip- 
ticity and PA. The time delay measurement and the lens mod- 
eling, combined with integral field spectroscopy of the lens in 
SDSS J0924+0219 might therefore prove extremely useful to 
map the mass-to-light ratio in the lens, by comparing the lens- 
ing and dynamical masses, to the light distribution infered from 
the HST images. 
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